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Abstract 

Let A be an n x n matrix with entries djj in the field C. Consider the following 
two involutive operations on such matrices: the matrix inversion I: A i— > *4 _1 
. and the element-by-element (or Hadamard) inversion J: aij i— ► a^ 1 . We study 

| the algebraic dynamical system generated by iterations of the product J o I. In 

the case n = 3, we give the full explicit solution for this system in terms of the 
initial matrix A. In the case n = 4, we provide an explicit ansatz in terms of 
CO ', theta-functions which is full in the sense that it works for a Zariski open set of 

initial matrices. Lhis ansatz also generalizes for higher n where it gives partial 
q ■ solutions. 

•l-H 

"3 ! 1 Introduction 

> 

A composition of two noncommuting involutions acting on square matrices generates 
sometimes an interesting dynamical system (more specifically: an algebraic dynami- 
cal system with discrete time). A good example can be found in paper [1] (see also 

references there), where block matrices were considered, A, B, C and D 

being themselves matrices n x n. The first involution consisted in taking the usual 

matrix inverse (of the whole block matrix), and the second one was the following block 

(A B\ (A C\ 
transposing: I „ n ] i— > ( _ ) . Such a system (at least when considered to 



x C D J \B D / 
within some natural gauge freedom) was shown in [1] to be a typical solitonic system 
solved by usual algebraic-geometrical methods. 

There is a natural desire to find solutions to more general (than just solitonic) 
dynamical systems. Note that the dynamical system generated by the following trans- 
formation: z i — , where z belongs to the unit circle in the complex plane, is obviously 
solvable but exhibits (as much obviously) a chaotic, not solitonic, behavior. 



I 



The aim of this paper is to present solutions to the following algebraic dynamical 
system (see [2] about its origin and some results for its particular cases). Let A be 
a n x n matrix whose entries belong to the complex field C. We consider two 
involutive operations on such matrices: the matrix inversion I: A i— > A^ 1 and the 
element-by-element (or Hadamard) inversion J: a^ \— > a" 1 . Our dynamical system is 
generated by iterations of the product J o I. 

The case n = 2 is trivial but serves us as a useful warm-up exercise. For the case 
n = 3, we give the full explicit solution for this system in terms of the initial matrix A. 
For the case n = 4, we provide an explicit ansatz in terms of one-dimensional theta- 
functions. This ansatz provides a full solution in the sense that it encompasses a Zariski 
open set of matrices. For the case n > 5, the same ansatz also works but encompasses 
only a subvariety of matrices of nonzero codimension. 



2 2x2 matrices 



We denote the matrix and its elements as A 



and the determinant of A 



Despite the triviality of this system, its solution supplies us with some tool that will 
play a key role also in the case n = 3. This is the multiplicative basis of evolution. 

a b 
c d 

as A = ad — be. The crucial point is that both transformations / and J of the following 
six values: 

a, b, c, d, A and (— 1) (1) 
are described in purely multiplicative terms. Namely, / act like this: 



a ^ dA _1 , b ^ (-1)6A _1 , c 
and J like this: 



(-l)cA- 1 , d i-> aA-\ A i-> A" 1 , (-1) i-> (-1); 



(2) 



a ^ a" 1 , b ^ ft" 1 , c ^ c~\ d ^ cT 1 , A ^ (-l)a" 1 6" 1 c" 1 d" 1 A, (-1) ^ (-1). (3) 

One can say that there are two matrices Tj, Tj e GL(6, Z) which act multiplicatively 
(in the following sense: the number 2 acts on a variable x multiplicatively by x i— > rc 2 ) 
on columns of values (1), where 



and Tj 



Of course, the answer for values (1) after iV steps of evolution is given in the similar 
way by the matrix (TjTi) n . This latter is given by two slightly different expressions 
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for N odd and even. For example, if N — 2k, then 



(TjTrf = 
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This means that the new value of a after 2k steps is 

a(2k) = a- k+1 b- k c- k d~ k A 2k (-l)- k } (4) 

and similarly for b, c and d. In the r.h.s. of (4), of course, the initial values of the 
variables are taken. 



3 3x3 matrices: one special formula 

This section provides a formula necessary for building the multiplicative basis for the 
evolution of 3 x 3 matrices, in analogy with section 2. 
For a matrix 

(a b c 
f 9 h\ (5) 
r s t 

we denote dh^4 the determinant of its Hadamard inverse multiplied, for convenience, 
by the product of all elements of A: 

dhA = agtbhr + fscagt + bhrfsc — rgcfbt — fbtsha — shargc. (6) 

It turns out that dh A behaves very nicely under the (usual) inversion of A: 

dhA' 1 = -dhA{detA)-\ (7) 

The direct proof of formula (7) consists simply in applying computer algebra. This, 
however, does not explain how one can arrive at such formula. In the rest of this section 
we present heuristic argument which clearly shows that a formula of such kind must 
exist. Our argument was suggested by paper [3]. 

Let dh^4 = 0. This means that the matrix 

fa- 1 b- 1 c- 1 
B = J{A) = f- 1 g- 1 h- 1 
yr" 1 s- 1 r 1 

is degenerate (remember that we do not care about the rigor!). This means that 

a\ I a 

B ( P I =0 for some (nonzero) column I (3 | . In terms of initial matrix A, this 

II \7. 
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yields: 

.A I (3' I A r = I (3" I . (8) 



fa' 











P' 
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Here the superscript T means matrix transposing; the values a', . . . , 7" are given by 

a R> P rJ 7 

a/r ogs cftt 

a" = a 'a 2 + (3'b 2 + 7 'c 2 , (3" = a' f + (3'g 2 + 7 'ft 2 , 7" = «V 2 + /?'s 2 + 7 't 2 . 

It is clear that a relation similar to (8) holds for 1(A) = A^ 1 as well. This means 
that (J o I) (A) = J (A' 1 ) is degenerate. In other words, 

dhA = dh(A)" 1 = 0. 

There is no need to make this argument rigorous because, as has been said, for- 
mula (7) admits a direct verification. 



4 3x3 matrices: the solution 



We proceed along the same lines as in section 2, which is possible due to formula (7) 
from section 3. We use notation (5) for the entries of matrix A. 

The multiplicative basis of evolution comprises now 21 values: the matrix elements 
a, b, c, /, g, ft, r, s, t; the determinant A = det A; nine cofactors of matrix elements 



denoted by corresponding capital letters: A 

the value S = dhA defined by (6) and, final! 
involutions act on these values. 
The matrix inverse /: 
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, B = 
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value (— 1). 


Here 


is 



and so on; 



a h 

A 

A 



AA' 1 , b h-> FA -1 , . . . , ftn HA' 1 ; 
> aA~\ B ^ /A -1 , . . . , if h-> ftA -1 ; 
>A-\ H^(-1)HA- 4 , (-!)-> (-1). 

The element-by-element inverse J: 

a 1— > a -1 , . . . , ft, 1— > ft -1 ; 

Ah> (-l)A# _1 ft _1 s _1 £ _1 , H i-> (— l)i?a -1 6 -1 / -1 ^ -1 ; 



A ^ Sa-^-^-V-^-^-V"^- 1 *- 1 , 
S i-> Aa-H^c^f^g^h-^s-H- 1 , 

(-1) -(-I)- 



A matrix TjTj corresponding to a step of evolution can now be calculated in analogy 
with section 2, but now it has sizes 21 x 21, and we do not write it out here. Still, 
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it makes no difficulty for a computer to handle such matrices. The remarkable fact is 
that all the eigenvalues of TjTj are sixth roots of unity, as the following table shows: 



eigenvalues 


l+v^3 
2 


1— /=3 
2 


2 


-1— /=3 
2 


1 


-1 


multiplicity 


1 


1 


4 


4 


7 


4 



Of course, XjT/ does have nontrivial Jordan boxes. An interesting thing with them is 
that they all correspond only to eigenvalues ±1. 

Here is the explicit answer for the variable a after 3/c steps: 

a(3k) = a -2fc+i fe -fe c -fej-V- fc A 2fe A- 2fe S- fc C- fc F- fe J R- fc S 2fc (-l)- 3fe ^- 1 ). (9) 

Again, in the r.h.s. of formula (9) the initial values of all variables are taken. 



5 4x4 matrices: the conservation laws and com- 
puter algebra results 

In this section we describe the way that has actually led to the ansatz presented in 
subsequent sections. Formally, however, the construction of our ansatz does not rely 
on the material of this section. 

First, the evolution of a 4 x 4 matrix A can be considered to within the following 
gauge freedom: we can consider not A itself but its equivalence class with respect to 
its multiplication both from the left and from the right by some diagonal matrices B 
and C. Clearly, if 

A' = BAG (10) 

then 

(J o I)(A') = C(J o I)(A)B. (11) 
Using a transformation (10), we can reduce (almost any) A to the following form: 

/I 1 1 1\ 

a b 

f 9 

r s 



A 



1 
1 

V 1 



c 
h 



(12) 



So, in a sense, there remains in A nine parameters. If we find eight conservation laws 
this will be a strong argument suggesting that the evolution goes just along elliptic 
curves (because only elliptic and rational curves have infinite number of automor- 
phisms) . 

Second, there are many conserved quantities which are, moreover, invariant un- 
der (10). Consider a decomposition of a 4 x 4 matrix in four 2x2 submatrices, for 
instance, one of the following: 

/o o * * \ / §°§\ 



o 
o 



or 



o 



o 
* 
o 
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where the entries denoted by the same symbol belong to the same submatrix. 

For any such decomposition p, we construct the value U p — the product of four 
corresponding minors. One can verify that under one step of evolution, the ratio of 
any two such values H p /H q goes into 11^/11^, where IT^ is the product of the cofactors 
of the minors of decomposition p (in the new matrix A). So, U p /Il q is invariant under 
two steps of evolution. 

Computer experiments show that there are eight algebraically independent invari- 
ants of such kind. Moreover, given fixed values of these invariants and using the 
form (12) for the matrix, one can exclude all except two parameters in (12) from the 
equations and get the curve as given by just one equation in two variables. Its genus 
turns out to be one, as expected. Besides, one can see from that equation that, for 
example, there are four points in the curve where the function a(z) (z being a param- 
eter for the elliptic curve) takes value 1, and the function b(z) takes value 1 in exactly 
the same points, as well as some information about the coincidence of some poles and 
zeros of those functions. 

This was exactly what led us to the ansatz presented in the following sections. 



6 A determinant of theta-function ratios 



K = (kij), where k %J = r^—^ - ( 13 ) 



The key formula for our ansatz is the formula for the determinant of the n x n matrix 

#(y - A, - fij) 
■d(x + \ + fi 3 ) ' 

Here % stays of course for the number of a row and j for the number of a column. So, 
there are complex variables x and y and, moreover, two arrays, (A«) and (fij), each of 
n complex variables. By -d we denote here the odd Jacobian theta-function: 



oo 

d{u) = 2q 1 ' A sin ^ ff(l - 2q 2n cos ^ + q 4n ){l - q 2n ), 
2K - LJ - K 

n=l 

where q = exp(—TcK'/K); K and K' are called half -periods. Our formula states that 

1 n n 

det/C = -($(x + y)Y $((n - l)x - y + + ^ pi^) 

i=l j=\ 

Yl ^(A i2 - AiJ • Yl - Mil) 

l<il<i 2 <« l<ji<?2<« (-\A\ 

n n ' V / 

IJIJtf(x + A i + ^) 

i=l j=l 

Proof of formula (14) goes by induction on n. Consider det/C as a function of A n . 
It is obliged to have the following form: 

n 7= i + \ n + ^j) 
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where F and G are some quantities that do not depend on A n but can depend on other 
variables. Here are the reasons why exactly such factors appear in formula (15), typical 
for proofs of the like formulas for theta-functions. 

Each factor i?(A n — Aj) in the numerator is responsible for the fact that det /C 
obviously vanishes when A n coincides with any other Aj (because of two identical rows). 
The denominator of (15) is simply the common denominator of all elements of /C. 
Finally, the factor i?(A n — G) is necessary to ensure that the whole expression has the 
same number of zeros and poles as a function of A n . 

On the other hand, in the neighborhood of the value A n = — x — /i n , where k nn has 
a pole, our determinant behaves as 

det JC « det /C smaUer H x + v) — ^ 

V(X + A n + fj^) 

where /C sma n er is the same matrix K, but without its nth row and nth column. Comparing 
(16) with (15) (where we at this moment also substitute A n = —x — /i n ), we can get 
the quantity F in the following form: 



&(x + y) n"=i + /ijQ 



F = det /C smallcr - „ ^n-i \ — ■ ( 17 ) 



Now, substituting (17) in (15), we get the following, almost final, formula expressing 
det K, through det K, STaa \\ er : 



det /C = det /C sma n cr 

+ y) Y\%1 0(-l*n + to) ^(A n - G) Y[ n ~l ^(A n - A,) 



(18) 



V{-X — fl n — G) UiZl -»n- A,) n;=l + An + to) 

It remains to calculate the value G. It can be deduced from the following reasoning. 
According to the inductive hypothesis, det/C sma n er contains the multiplier -#^(n — 2)x — 

V + YllZi Aj + Sj=i to^j ' so ^ mus ^ have a zero at such A n when the argument of that 
theta-f unction equals zero. On the other hand, det /C, generally, does not have a zero 
at such A n . So, the mentioned theta-function must cancel with the same factor in the 
denominator of (18), which role only ■&(— x — ji n — G) can assume. This implies 



n-l 



-x - fi n - G = ±((n - 2)x - y + ^ Aj + ^ n^j , 

i=l j=l 

and the sign here, namely plus, can be fixed for example by considering the expres- 
sion (18) as a function of x (which must have the right difference, namely n, between 
the number of poles and zeros in a parallelogram of periods). 

Once we have got the proper formula for transition from det /C to det /C sma n er , the 
inductive step is over. As for the induction basis, the formula (14) does obviously hold 
for n = 1. 
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7 Expression for a matrix element of ( J o I) A 



Let us apply the transformation Jo/, i.e. one step of our evolution, to the nxn matrix 
A = K, given by formula (13). A matrix element a^ ew of the obtained matrix ^4 ncw is 
given by formula 

aT = (19) 

where Aij is the cofactor for the element of A. Both the numerator and denominator 
in (19) are determinants of the form (14). The calculation yields: 

new global row column clement (Of]} 

a ji — u u j u i u ij i l zu J 

where 

a 



n n 

global 



= $(x + y)$((n-l)x-y + J2 X * + J2l 2 i) ( 21 ) 



k=l 1=1 

is a factor which depends neither on j nor on i; 



nLi + x k+ Hi) 



a T W = A" ,7 ^ ( 22 ) 



depends only on j (the row number for a°f w ); 



j 

n—l 



, column 



nr=i^+^+^) 

depends only on i (the column number); and the last factor 



(23) 



^element = + K + N) ^ 

$((n -2)x-y + Y2=i X k + Er=i Vj) 

has much the same form as the initial a^. 

To be exact, the difference between the matrix (a°' emcnt ) made of matrix elements 
(24) and the initial matrix A = /C can be described as follows: change 

n n 

i=i j=i 

and then perform the matrix transposing. As for the factors (22) and (23), their effect 
consists in multiplying the matrix (a°' emcnt ) from two sides by diagonal matrices, i.e. 
doing a gauge transformation (10). The main point is that if we have done not one 
but N steps of evolution, the effect of all arising factors (22) and (23), as well as (21), 
consists just in the appearing of some products of theta-functions with their arguments 
changing according to a simple law. We do not write out here the corresponding 
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obvious but bulky formulas. What we see already is that the evolution of a matrix of 
the form (13) can be described by an explicit formula. The same applies, obviously, to 
a matrix that was obtained from such one by a transformation (10). 

Remark. Here we do not mean just the evolution of gauge equivalence classes as 
in section 5. Formula (11) shows that if we know what happens with A, we also know 
what happens with BAC after any number of evolution steps. 

8 Comparing different values of n 

8.1 n = 3 

The ansatz (13) (taken together with the possibility of multiplying a matrix by two 
diagonal matrices as in formula (10)) is definitely superfluous for the case n = 3 where 
we have presented, in sections 3 and 4, a more direct approach which, by the way, 
gives the exhaustive information about the cases where the evolution cannot go ahead 
because of a division by zero. 

8.2 n = 4 

Here the ansatz (13) together with multiplication by two diagonal matrices contains 
exactly 16 independent parameters, i.e. gives a Zariski open set of matrices. To see 
this, note first that the 10 values x, y, Aj and /ij produce really only 8 parameters, 
because nothing in (13) changes if we do one of the following translations (a being an 
arbitrary complex number): 

Aj i— > Aj + a for all i, /ij i— > fij — a for all j 

or 

x >-> x + a, y i— > y — a, Aj i— > Aj — a for alii 

Second, the modulus of the elliptic curve is the 9th parameter. And finally, the two 
diagonal matrices produce the 7 remaining parameters. 

8.3 n > 5 

Our ansatz works for any n, but when n > 4 it corresponds only to a subvariety of 
a nonzero codimension in the space of all n x n matrices. For instance, for n — 5, 
the calculation of parameters similar to that done in subsection 8.2 shows that we can 
describe in the same way the evolution of a 20-parameter family of matrices. 

9 A special case with period 4 

Although the parameterisation (13) together with the possibility of multiplying the 
matrices by two diagonal ones as in (10) encompasses in the case n = 4 a Zariski open 
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subset of all 4 x 4 matrices, it does not include some interesting special cases (and it 
does not seem very easy to obtain them as any limiting cases). One specific feature of 
parameterisation (13) is that if some 2x2 minor of the matrix equals zero, then all 
other minors also do so. We present here a matrix 



A = 



(I 
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1 
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xy 
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zt 
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V 1 


t 


y 





(25) 



that does not obey such a requirement. Its interesting property is that if we consider 
its evolution up to gauge transformations (10), then after four steps we get back at the 
initial matrix (25). This statement is proved by a direct computer calculation. 
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